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A self-consistent theory of the frequency dependent diffusion coefficient for the Anderson 
localization problem is presented within the tight-binding model of non-interacting electrons 
on a lattice with randomly distributed on-site energy levels. The theory uses a diagrammatic 
expansion in terms of (extended) Bloch states and is found to be equivalent to the expansion 
in terms of (localized) Wannier states which was derived earlier by Kroha, Kopp and Wolfle. 
No adjustable parameters enter the theory. The localization length is calculated in 1, 2 
and 3 dimensions as well as the frequency dependent conductivity and the phase diagram 
of localization in 3 dimensions for various types of disorder distributions. The validity of a 
universal scaling function of the length dependent conductance derived from this theory is 
discussed in the strong coupling region. Quantitative agreement with results from numerical 
diagonalization of finite systems demonstrates that the self-consistent treatment of cooperon 
contributions is sufficient to explain the phase diagram of localization and suggests that the 
system may be well described by a one-parameter scaling theory in certain regions of the 
phase diagram, if one is not too close to the transition point. 


1. Introduction 

The concept of localization of a quantum particle in a random potential was 
introduced by P.W. Anderson [1] in his pioneering work in 1958. Although 
considerable progress has been made since then, the problem still resists a complete 
theoretical understanding. In further developing ideas by Thouless [2], Abrahams 
et al. [3] formulated a real space scaling theory of localization in 1979, which is 
based on the assumption that the dimensionless conductance g(L ) of a sample with 
finite length L is the only relevant scaling parameter. The assumptions of the one- 
parameter scaling theory were supported by the work of Wegner [4] and others 
[5-7], who mapped Anderson’s original Hamiltonian onto non-linear a models of 
interacting matrices and applied perturbative renormalization group calculations in 
d = 2 + e dimensions. Self- consistent theories [8, 9], developed at the same time, 
also yielded scaling [10] in agreement with these field theoretical treatments. In this 
way it was established that the Anderson transition is a continuous phase transition 
characterized by a diverging length, critical exponents and an order parameter. 


However, in recent years the validity of the one-parameter scaling theory has 
been called into question by several new developments: (1) The discovery of large, 
universal conductance fluctuations in mesoscopic samples [11, 12] has raised the 
question whether the distribution of conductances rather than the ensemble-averaged 
conductance g(L ) would obey one-parameter scaling [13, 14] or whether more than 
one scaling parameter would be required in the strong coupling regime [15, 16]. 
(2) Recent results [17, 18] concerning localization on a Bethe lattice, which were 
obtained within the supersymmetric matrix model introduced by Efetov [19], show 
a non-power law singularity at the transition point with an exponentially vanishing 
diffusion coefficient at the metallic side and a diverging localization length with 
critical exponent v = 1 at the insulating side. (3) According to Kravtsov, Lerner 
and Yudson [20], in the non-linear a model there exist relevant operators from higher 
gradients, which may indicate a violation of one-parameter scaling. Furthermore, 
Wegner [21] performed an e expansion for the non-linear a model and found in 
four-loop order a correction to the critical exponent v in d = 3 dimensions which 
violates an exact relation for v established by Chayes et al. [22]. Thus, it is unclear 
whether the asymptotic e expansion about d = 2 can be extended to 3 dimension^]. 
In this spirit, Zirnbauer [23] has done a non-perturbative calculation in d = 3 within 
Efetov’s model using the Migdal-Kadanoff approximate renormalization scheme. 

The question of the precise critical behavior near the Anderson transition remains 
controversial up to now. In the present paper we want to address a different problem, 
namely how to connect the critical regime with the parame- ters of the Hamiltonian 
in a quantitative way. To this end, the self-consistent theory originally devised by 
Vollhardt and Wolfle [8] is extended to the lattice and fully renormalized in the 
sense that the calculation of all quantities entering the self-consistency equation is 
extended to the strong coupling regime. As will be seen below, the self-consistent 
expansions in terms of small and large disorder, respectively, are equivalent at 
the level of maximally crossed diagrams. One-parameter scaling is an inescapable 
consequence of the diagrammatic expansion, unless there exist as yet undiscovered 
infrared-divergent contributions besides that class of diagrams. However, by 
comparing with independent numerical renormalization group calculations [24] for 
the same model one may gain some insight in how wide the critical region around 
the transition point is in which deviations from one-parameter scaling occur and in 
which other infrared singularities not contained in the present theory might become 
important. 

In the next section, quantities which behave non-critical at the Anderson 
transition are calculated in single-site approximation (CPA). Section 3 contains the 
solution of the Bethe-Salpeter equation for the density correlation function following 


1 See the discussion in ref. [21] 


Vollhardt and Wolfle [8], which leads to a frequency dependent diffusion coefficient 
D{uj) expressed in terms of the irreducible particle-hole vertex part U P)P '. In section 
4, U P)P ' is calculated by classifying all irreducible particle-hole diagrams in terms of 
the most important infrared-divergent contributions, and a self-consistent equation 
is derived in this way. From this equation, the f5 function of the length dependent 
conductance g(L) is obtained in section 5 along with a discussion of its validity in 
the strong coupling regime. Some concluding remarks are given in section 6. 

2. Calculation of non-critical quantities 

The model under consideration is defined by the Hamiltonian 

H = y^t n m\n){m\+y^y n \n){ri\ , ( 1 ) 

m,n n 

with random site energies V n , distributed according to probability P(V n ) and 
hopping amplitudes t nm from site m to site n [t nn = 0). Later evaluations will 
be done for isotropic nearest neighbor hopping t^ nm ^ = t and for box, Gaussian and 
Lorentzian probability distributions Pg, Pq and Pg, respectively, 

±Q(W/2-\V\) , 

75P7i exp (^) ' (2) 

W 

7r(V 2 + W 2 ) ' 

Here 0 is the step function, and the width W plays the role of the disorder parameter. 
The width of the Gaussian distribution is scaled in such a way that the second 
moments of Pb and Pq coincide. 

The usual expansion with respect to small disorder W/t is done in terms of Bloch 
states repeatedly scattering off the impurity levels (extended state expansion). Thus, 
the single-particle Green function is expanded as 

G^ a (E) = (OKE-HPiO)- 1 ^) 

= G^ (0) + E G ^ (0) ^ (0) 0) 

+ E G n/ A(0 V i Gg M(0 V J Gf T / A(0) + ... , 


Pb(V) = 
Pg(Y) = 
Pl{V) = 






where 


1 


qR/ a (°) 


(E) 


E — e p ± iO 


(4) 


is the retarded/advanced free Green function of a Bloch state in momentum 
representation and e p = 2t Yli=i cos pta is the dispersion on a d-dimensional simple 
cubic lattice with nearest neighbor hopping. 

When the Green function G, averaged over the ensemble of all realizations of the 
random potential, is calculated, 


G r ' a (E) = ( G r ' a {E)) = [il [dV n P(y n )\ G r ' a (E) , ( 5 ) 

n 


impurity levels V which happen to be at the same site have to be averaged coherently. 
This is expressed diagrammatically (fig. 1) by attaching the corresponding impurity 
lines to the same point and introduces irreducible diagrammatic parts in the 
calculation of averaged quantities. Then, G satisfies the Dyson equation 

G = G^ + G^EG , ( 6 ) 


which is easily solved in momentum space: 


G* /A (E) = 


E- £p -Z* /A (E) 


(7) 
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Figure 1. Construction of MOCs: dashed lines represent interactions with the impurity 
potentials V, horizontal double lines the full Green function G. Each MOC generated from 
an irreducible diagram (left) corrects for coinciding sites in the respective ’’direct” diagram 
with the same topology. 











The self-energy S p '‘ (E) represents the sum of all irreducible diagram parts with 
respect to cutting one propagator line G . The averaging process causes a technical 
difficulty which is due to the discreteness of the lattice and, therefore, is present in 
any lattice model: In any averaged diagram, lattice sums are restricted to sites 
which do not occur anywhere else in the diagram, since coinciding sites would 
change its topological structure. This difficulty is lifted exactly by extending the 
sums over all lattice sites and correcting for the oversummations by subtracting 
”multiple occupancy corrections” (MOCs) [25]. MOCs are generated out of any 
given irreducible diagram by breaking off impurity lines from common sites in all 
possible ways and subtracting these terms from the original diagram. MOC diagrams 
may, in turn, require MOCs (next generation) if they have multiple impurity lines 
at any site. The construction of MOCs is illustrated in fig. 1. 

Alternatively, the Hamiltonian (1) allows an expansion in powers of t/W, treating 
the kinetic term as a perturbation (locator expansion) [25, 26]. The unperturbed 
on-site Green function, or locator at site n, is then given by 


9nL A (E) 


E — V n ± iO 


( 8 ) 


Diagrammatically, each term in this expansion is represented by (vertical) locator 
lines g attached to the corresponding sites and connected by (horizontal) hopping 
lines t. Similarly as in the extended state expansion, it is useful to define a self¬ 
resolvent S(E) as the sum of all irreducible diagrams with respect to cutting one 
hopping line, so that the full, averaged Green function satisfies the Dyson equation 


or, in momentum space, 


G = S + StG , 


G« /A (E) 


1 



( 9 ) 

( 10 ) 


By comparison of eqs. (7) and (10), we can now establish an important 
relation between the irreducible parts of the weak and strong coupling expansions, 
respectively [27]: 


£p /A {E) =E - [Sp^iE)}- 1 . (11) 

Let us now turn to the actual evaluation of single-particle quantities. This will 
be done in single-site or coherent potential approximation (CPA), which amounts 
to summing up all irreducible diagrams, including MOCs, that do not contain any 




crossings of lines 
by 


where the brackets denote the ensemble averaging and Go = (1/iV)^ G p . Note 
that in CPA the self-energy is momentum independent. The CPA can be viewed as 
the first order of a self-consistent cluster expansion and is known [28] to become 
exact in the limit of large lattice coordination numbers (or high dimensions). 
Furthermore, by employing the identity (11) it can be shown [29, 27] that the 
single-site approximations for the extended state and for the locator expansion are 
equivalent. Thus, the CPA interpolates between the exact limits of strong and 
weak disorder. (It cannot account for a correct description of the singularities at 
the band edges, however.) Therefore, since averaged single-particle quantities are 
smoothly varying functions at the Anderson transition [30], the CPA is expected to 
be sufficient for calculating these quantities (fig. 2). 

Kopp [31] formulated a version of the Vollhardt-Wolfle theory [8] using the locator 
expansion. While, however, this formulaton was not invariant with respect to a shift 
of the energy scale, a fully invariant formulation is given in ref. [27]. It will be shown 
below, that the latter coincides exactly with the formulation using the extended state 
expansion. Therefore, we will confine ourselves to the extended state formalism in 
the following. 


. The resulting self-consistency equation for the self-energy is given 


V - T, 


R/A 


(E) 


1- [V -Zo /A (E)]Go /A (E) 


= 0 


( 12 ) 



Figure 2. Disordered density of states in CPA for Gaussian distribution. Data points 
taken from ref. [24] for comparison. 




3. Correlation functions 


In order to investigate the localized or extended character of states, four- 
point functions have to be considered, since averaged single-particle functions 
are translationally invariant. The averaged particle-hole (p-h) Green function 
<Ppp,(E,oj,q) = (G r t (E + uj)G a ; (E)}, with p± = p ± q/2 etc., is calculated 

from the Bethe-Salpeter (B-S) equation 

&${E,u>,q) = GK + (E + u>)G$_(E) 5 pp , (13) 

+ G$ + (E + u)G£jE)±^U$,(E,u,q)^(E, U ,q) , 

p" 

where denotes the p-h irreducible vertex function. The B-S equation is 

supplemented by an exact Ward identity, valid for all frequencies [8]: 

(E + u) - St = U p^ G p + (E + u)- G$_ (E)} (14) 

p' 

Before we will attempt to solve the B-S equation, it is useful to consider the 
density correlation function <P pp = (l/N 2 ) ]T) pp , i n the simplest, i.e. ladder 

approximation. That is, the irreducible vertex U is calculated as the sum of all 
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Figure 3. Diagrammatic definition of the single-site vertex U^ A , including MOCs. 





















irreducible p-h diagrams without crossings (single-site approxmation, fig. 3). Note 
that this class of diagrams is generated from the CPA diagrams of section 2 via the 
Ward identity. The result may be obtained by direct summation (in analogy to the 
procedure shown in ref. [27]) or by making use of the single-site version of the Ward 
identity, E* - ^ = U RA (G R - G A ), 

RA V a (E)-^{E + u) 

° G$(E)-G*(E + u>) ( } 

The ladder approximation for the total (reducible) vertex is then easily evaluated 
as (fig. 4) 


F L(u,q) = 


ttRA 

u o 


o,q ->o 2-KlN(E)(U q A ) 2 


1 - U RA (l/N) E P G r + G a _ uj + iq 2 D 0 

where N(E) = (l/7r)ImG^ is the density of states and 

1 1 


(16) 


Hn = 


ttN(E) N 


X> P -q) 2 ( Im ^ 


A\2 


(17) 


is the ’’bare” diffusion coefficient (v p = Ve p ). For the evaluation of T^ in the 
hydrodynamic limit, uj, q —>■ 0, we have used the expansion 


1 _ s~iR 

G r + {E + u)G a _{E) = 


X 


0 

1 - 


- §(Gp ~ G p)( v P -q)9 ( 18 ) 


U ra (G a -G r ) 


1 g a -g r 

In ladder approximation, the density correlation function reads 


c p L = 


pp N 




oj , q ^o 2mN(E) 


co + iq 2 D 0 


(19) 


The expected diffusion pole structure of (19) is a consequence of the Ward identity 
(14) and indicates particle number conservation. The diffusion coefficient D is 
expected to vanish at the Anderson transition. A perturbational calculation for 
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Figure 4. Ladder approximation for the total p-h vertex. A full vertical line represents 
the single-site vertex Uq. 


D, using, e.g., the standard Kubo formula, would therefore require to take into 
account arbitrarily small contributions, i.e. infinite orders of perturbation theory. 
This formidable task can, however, be circumvented by observing that a vanishing D 
is equivalent to a divergence of the density correlation function <P pp at ui —>• 0, q / 0. 
Therefore, for uj —>• 0 it is more useful to identify and take into account the largest 
infrared divergent contributions to ( P pp , which will enter through the irreducible 
vertex 17(^1 • This will be done in section 4. Here, we first solve the B-S equation in 
the hydrodynamic limit generally in terms of 

Eq. (13) can be rewritten in form of a kinetic equation by using 


G 


R 

P+ 



r<A _ r^R 
Lt p- ^ 7 p+ 

1 _1 

QR qZ— 

u p+ °p- 


multiplying the denominator to the left-hand side and summing over p', 


1 1 


" - (v P ■ q)9 + = K ~ °?)( 1+ -rn E < 20 > 
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N 2 


Using again the Ward identity (14) and summing over p' explicitly incorporates the 
conservaton of particle number in the form of the continuity equation 


u<P pp - q<P Jf) = 2iriN(E) . (21) 

It relates <P pp and ( P Jp = (1/N 2 ) ^ pp /(v p • q)^^(q), the density-density and 
the current-density correlation functions, respectively, to each other. It is this 
conservation law that governs the hydrodynamic behavior of the system. Therefore, 
<P p = (1/N) <Ppp, may be approximated by its projection onto density and 
current correlations only [31], 


<Pp = Ap<P pp + B p (q)& jp + A<P p . 


( 22 ) 



















Eq. (22) can be viewed as an expansion by moments of the current vertex (v p • q). 
In fact, contributions A^ p from higher moments are less divergent and may be 
neglected in the hydrodynamic limit. Since the critical properties of the localization 
transition are contained in the relaxation functions (£ pp and <Pj p , the coefficients 
A and B behave uncritical and may now be obtained from the simple ladder 
approximation, where <P p , <P pp and <Pj p are given explicitly, 

IrnG'p 

(l/A0£ p ImC£ ’ 

(v p • q)(ImG^) 2 

(IfflEpIvp'tfllniCJ) 5 ' p3) 

We use eq. (22) to decouple the momentum integrals in the kinetic equation (20) 
and obtain, after multiplying with (v p -q) and summing over p, a current relaxation 
equation [8], 


^4p — 
^p(q) = 


co + 2iImE^ + iM(cu) <Pj p — q ■ 2ImS^ Do <P pp = — qR(E ) (24) 

where the ”current relaxation kernel” M(co) is given by 

q) Im ^p t^pp' (imGp/) 2 (v p / • q) (25) 


— iM(co) = 


ttN(E)D 0 N 2 


pp' 


and 


R(E) = ^ E( v p-q) 2 ^ G p 2 + G : 

p 


R2\ 

P > 


(26) 


The density correlation function follows from (21) and (24), 

2«N(£) - g G-HllMMMI 


®PP = 


(27) 


co + iq 2 D(co) 

Here, we have introduced the generalized, frequency dependent diffusion coefficient 

Do 


D(co) = 


1 — icoT + tM{co) 


(28) 


where 1/r = 2ImS^ is the single-particle relaxation rate. The ladder diffusion 
coefficient Do has been renormalized by a factor essentially given by the inverse 










current relaxation kernel, which is in turn expressed in terms of the irreducible 
vertex [7^1. Thus, in the hydrodynamic limit it is indeed sufficient to take into 
account the strongest infrared divergencies in , and they will eventually make 
D(uj —>• 0) vanish at the transition point. Also because of this divergence, the term 
~ q 2 R(E ) in eq. (27) may be dropped for u>, q —>• 0. 


4. Infrared divergence of the irreducible vertex and self-consistent 
theory 

It was shown [8] that the infrared divergence of the diffusion ladder T^ as an internal 
part of the total irreducible vertex U^, is cancelled by vertex corrections There is, 
however, another class of diagrams that introduces such a singularity into U^,. 
namely the class of maximally crossed diagrams known from ’’weak localization” 
[32, 33] (fig. 5a). Unlike the diffusion pole, this singularity has no classical analog. 
It is only present in systems with time reversal symmetry. Due to this symmetry, the 
p = — p' singularity (”2 kp singularity”, ’’cooperon”) present in the particle-particle 
channel is carried over to the p-h channel. In this way, the sum of maximally crossed 
diagrams is obtained as 


A pp ' = T L (Q) 


2t TiN(E){U* A ) 2 
w + z(p + p') 2 A) 


(29) 


Here we have introduced Q = p + p'. 

In any p-h vertex diagram containing a A block the singularity at p = — p' 
will be strongest if A is crossing the total diagram diagonally, since otherwise the 
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Figure 5. (a) Sum of maximally crossed diagrams A pp /. (b) Exact classification of all 
diagrams of the irreducible vertex in terms of cooperons (wavy lines). 

























singularity is integrated over and thus weakened. Therefore - provided that there are 
no other singular contributions besides cooperons - the infrared divergent behavior 
of each diagram contained in the irreducible vertex ifejy can be classified by the 
number of diagonally crossing cooperons as shown in fig. 5b (Vollhardt and Wolfle 
[8, 33]). Wavy lines denote cooperon blocks, while the internal part is the sum 
of all reducible and irreducible p-h diagrams which do not contain any diagonally 
crossing interaction line. U represents all irreducible p-h vertex diagrams with the 
same restriction. 

By flipping over the hole line in the last diagram of fig. 5b and employing time 
reversal symmetry, the maximally crossing lines A pp / are disentangled into ladders 
Tl(Q) with the momentum argument q replaced by Q = p + p'. In this way it is 
seen that consists of all diagrams contributing to the density correlation function 
d> pp (uj,Q) except those with ladder diagrams at the ends: 


$(Q) = $ PP -R(Q)r L (Q)R(Q)-R(Q)r L (Q)3(Q) 

${Q) r l (Q) R(Q) - R(Q) r L (Q) $(Q) t l (q) r(Q) , ( 30 ) 

with R(g) = G p+ G p_ , 

p 

or 

= (Sw= (bfei) l ^ W) - * 2( « )riW)l ■ < 31 > 

From the last equality it can be seen that the divergencies of the crossing cooperons 
in the last diagram of fig. 5b are cancelled by the factor T l{Q) 2 in &(Q), leaving 
the singularity of ^ pp itself. Furthermore, diagrams with less than two diagonally 
crossing Cooper blocks A (third-fifth terms in fig. 5b) are less divergent and may 
be dropped. Thus we have from the second and the sixth diagram: 

= r L (Q) + T L (Q) <P(Q) T l (Q) = (U* A ) 2 d> pp (Q) . (32) 


^pp)’ governed by the pole structure of <P pp , which in turn feeds back into the 
diffusion coefficient eq. (28). In this way a self-consistent equation for D(u>) is 
obtained: 


D(uj)( 1 — icor) 


D ° + I]( v p • q) 


pp' 


IrnG’p (ImGp,) 2 
'—iw/D(uj)\ + (p + p') 2 


( v p' • q) 


2ImX$ 
[irN(E)] 2 D 0 ' 


with A 


(33) 






It is remarkable that this equation is identical to the one derived in the formulation 
in terms of a basis of Wannier states [27]: The expansions with respect to small 
disorder and with respect to small hopping amplitude are equivalent if the effects 
of quantum interference (cooperons) are taken into account in both formalisms in a 
self-consistent way. 

Because the cooperon propagator connects momenta of equal size and opposite 
direction, the term in eq. (33) renormalizing the bare diffusion constant is always 
negative. In dimensions d < 2 all states are localized even for arbitrarily small 
disorder, while for d > 2 one finds a continuous metal-insulator transition with a 
scaling relation [8] v = s/(d— 2) for the critical exponents s, v of the conductivity and 
the localization length, respectively, in agreement with field theoretical treatments 

[4-6]. 

Eq. (33) has been solved numerically, where the momentum integrals have been 
evaluated using the isotropic energy-momentum relation that leads to the model bare 
density of states Nq(E) = Nd(B/2)~ [ [1 — E 2 /(B/ 2) 2 ]” 1+d / 2 , i.e. that is defined by 
A r o(e p ) = [Sd/(2ir) d ]p d ~ l \dspfdp\~ 1 . Here B = 2Zt is the bare band width for 
lattice coordination number Z, Sd is the surface of the d-dimensional unit sphere 
and Ni = 1/n, IV 2 = 1/2, IV 3 = 2/ir. In d = 1 this approximation becomes exact. 

In fig. 6 the phase diagrams in the W-E plane are shown for box, Gaussian and 
Lorentzian level distributions P(V). The calculated phase boundaries are in very 
good agreement with results of exact numerical diagonalization [24]. In particular, 
the re-entrant behavior for energies outside the bare band is reproduced. This 
behavior is due to the fact that for energies near the band edge the density of states 
and thus the diffusion rate first increases with increasing disorder before localization 
effects start to dominate and drive the system to the insulating phase. The existence 
of two phase transitions at energies greater than B/2 = Zt can also be seen in fig. 7, 
where the u = 0 conductivity is plotted as a function of disorder for various values 
of the energy E. 

In fig. 8 the localization length £, defined by lim^_;.o[ — iu/D{u)] = l/£ 2 , is shown 
in the band center (E = 0) as a function of disorder in d = 1, 2, 3. The agreement 
with the numerical results [24] is again very good for d = 1, 2. For d = 3 deviations 
occur, which may be due to the fact that an accurate numerical determination of 
the critical disorder W c is crucial in this case. Asymptotic power laws in the regions 
of large and small localization length are also shown in the figure. In particular, 
from eq. (33) it follows that £ ~ IE -1 / 2 for IE —>• 00 . This is the expected result 
in a tight binding model because of the following argument: In the limit of small 
localization length of a wave function located at site i, £ will be determined only 
by the potential step of order IE at the nearest neighbor sites and will be equal to 
the decay length of the wave function into that potential step. For massive particles 
(which are always implied in a tight-binding model; mass m) this decay length is 







w/t 

Figure 7. DC conductivity in d = 3 dimensions normalized to the CPA conductivity 
c7o(uj = 0) as a function of disorder W for various energies E (box distribution). 


given by 1 /-y/2 m(W — E ), thus explaining the VF^ 1 / 2 behavior for large W. The 
same result is obtained from a diagonalization of the corresponding two-site problem 
in a Wannier basis. As seen from fig. 8, the crossover from the scaling behavior near 
the transition to the atomic limiting behavior indeed occurs at localization lengths 
smaller than one lattice constant a. Note that the numerical data show a similar 
crossover. 

In fig. 9 the real part of the dynamical conductivity in three dimensions is shown 
for box-shaped disorder distribution. In the metallic phase (W < W c ). Re<r(u;) has 
a dip at ui = 0 with a ~ y/Lo singularity, which eventually makes Recr(0) vanish at 
W = W c with a singularity ~ lu 1//3 . In the localized regime (IF > We c ), a low- 
frequency behavior Recr(u;) ~ ui 2 is obtained. For high frequencies cr(u) merges to 
the Drude behavior; the conductivity sum rule is fulfilled. 


<r- Figure 6. Phase diagram of localization in d = 3 for (a) box, (b) Gaussian, (c) Lorentzian 
disorder. M: metallic, I: insulating regions. Phase boundary according to eq. (33) (solid 
line) and ref. [24] (dots). The error bar of ref. [24] at E = 0 is indicated: ±0.5 1. Also 
shown are the band edge as determined in CPA (solid line in (a), (b)) and numerically 
in ref. [24] (short-dashed line) and the exact upper bound (long-dashed line in (a)). The 
critical disorder strength W c for localization in the band center (E = 0), as obtained from 
the self-consistent theory, is given for each case. 










Figure 8. Localization length in units of the lattice spacing as a function of disorder in 
dimensions d = 1, 2, 3 (a)-(c), for half band filling and box distribution, calculated from 
eq. (33) (solid lines) and according to ref. [24] (dots). Asymptotic power laws for weak and 
strong localization are indicated. 






Figure 9. Dynamical conductivity Re ct(lo) normalized to ero(O) versus normalized 
frequency lot for various values of disorder box distribution. Insets show the low- 
frequency behavior. 


5. Scaling function in the strong coupling region 

The controversy between recent results concerning the scaling behavior at the An¬ 
derson transition (compare section 1) and the fact that one-parameter scaling is 
observed in numerical renormalization group calculations [34, 35] at least for ener¬ 
gies close to the band center has led the author to re-examine the (3 function for the 
length dependent conductance derived from the theory presented above. The deriva¬ 
tion follows the discussion by Vollhardt and Wolfle [10]; however, special emphasis 
is put on the validity of the approximations made in the strong coupling region. 

Density response in a finite sample 

According to Thouless [2] and Abrahams et al. [3] the dimensionless conductance 
g(L ) of a sample with finite length L, g{L) = ( e/h ) 2 L d ~ 2 <r(L), should be a relevant 
scaling parameter. We want to calculate g(L) directly from a density response 
theory. 

Due to the broken translational invariance in a finite sample, the response is 
non-local in momentum space, so that the transport coefficient cr{L) cannot be 












calculated in a trivial way from the diffusion coefficient D(L) as in the infinite 
sample. Instead we consider the density response to an externally applied potential 
U(x) = e|E|(L — x) in position space: 

5p(x) = f dx'x{x — x',oj = 0) U{x') . (34) 

Jo 

Here x( x — x',u) is the one-dimensional Fourier transform of the density response 
function 

%(q, w) = ^ &pp(q, u>) + N(E) N(E) + l) , (35) 

with £ = lirn^^o [iD(uS) / , the localization length of the infinite sample. In this 

way, one obtains without approximations 

= e\E\N(E) hl + y)e~y , (36) 

x=L 1 

where y = L/£, the sample length scaled to the characteristic length £, has been 
introduced. This density gradient at the end of the sample gives rise to a diffusion 
current jdiff = — eD(L ) dp/dx\ x= L, which in the stationary case is equal and oppo¬ 
site to the electrical current j e i. 

Length-dependent diffusion coefficient 

The effect of a finite sample is that all wave numbers in the system are restricted 
to values greater than 1/L. Thus, the length-dependent diffusion coefficient D(L) 
can be obtained in analogy to eq. (33) as 

D(L) = D 0 + aT £ (VQ- P .q). (37) 

IQIM/F p /? V 

The main contribution to the p integral in (37) comes from the peak of the imaginary 
part of the Green function at the Fermi momentum, p = pf> Q — p = Pf- Also, as 
will be seen below, D(L) can be expressed as an integral over the diffusion pole with 
Q values restricted to Q < 1/L. Therefore, the p integration in (37) can safely be 
extended over the complete Brillouin zone, if pfL > 1. Furthermore, if the sample 
size is much larger than the mean free path, the localization length will be the same 
for the finite sample and for the infinite system. Thus, in eq. (37) the localization 


dp 

dx 





length £ of the infinite sample has been introduced. It is determined in the localized 
regime by eq. (33): 


0 = A) + Y ( v p ‘ q ) 


ImGp (ImGq_ p ) 2 


Q, p 


i/e 2 + q 2 


(vq-p • q) • 


Subtracting (38) from (37) we obtain 

D ( L ) = Y ' q ) 


IQI<i/i. p 


ImGp (ImGq_ p ) 
l/£ 2 + Q 2 


(vq-p • q) 


(38) 


(39) 


Now the two integrations factorize if Q < 1/L ^ l/£, where 1/7 is the inverse mean 
free path measuring the peak width of the Green function. After some algebra one 
obtains 


D{L) 


N(E) 


Do 7T \2nJ J o 


d y' 


,/ d — 1 


i + y 


r 2 ’ 


(40) 


where D$ is defined by 

= <j )2 ( im G )3 . (4i) 

u p 

Eqs. (40) and (36) together yielld the dimensionless conductance in the localized 
regime: 

t\ rl/y f i/d— 1 

g( L ) = c d y d ~ 2 { 1 + y) e~ y ^ dy' ^ +yl2 , (42) 

with the dimension-dependent constant c d = (l/pi)[Sd/{2 7r) rf ]. In the metallic 
regime (d > 2) one defines the characteristic length £ in analogy to the localization 
length [10], i.e. £ ~ (1 — X/X C rit) 1 ^ d ~ 2 \ and obtains 


9{L) = w 0 Ci T = 2 11+ r W 2 ) r < 2 -< i / 2 > y d ~ 2 \ 


(43) 


Clearly, the quantity g(L) = g(L)/(Do/Do ) depends only on the scaled sample 
length y and, therefore, obeys a universal scaling law under the conditions posed 
above. The corresponding /? function is explicitly calculated from eq. (42): 


P(g) 


ding =d _ 2 ( 1 + ^) e V _ V 2 

d In L g 1 + y 2 1 + y 


(44) 













Figure 10. Scaling function /3(g) (see text) in d = 1, 2, 3 dimensions. The dot indicates 
the zero of the /3 function in d = 3 as determined numerically in ref. [35]. 


where y = y(g ) from eq. (42). The function /3(g) is shown in fig. 10. 

Summarizing this section, we have found that a one-parameter scaling theory can 
still be derived from the self-consistent theory in the strong coupling regime. The 
validity of the one-parameter scaling assumption is however subject to the conditions 
PfL /$> 1 and L/£ 1, which were crucial in the above derivation. They set a 

lower limit for the sample sizes for which scaling can be expected to be observed. 
In particular, one-parameter scaling must break down at energies close to the band 
edge, where the Fermi momentum pp approaches 0. This is in agreement with results 
of numerical calculations [24, 36]. Furthermore, we find that only a ’’renormalized 
conductance” g(L) obeys scaling for large disorder. The renormalization term, i.e., 
Ing(L) — lng(L) = ln(.Do/-Do), which is the shift of (3(g) relative to (3(g), has been 
evaluated numerically. It was found that \w(Dq/Dq) = —0.10 at the transition point 
in the band center and that it varies only by a few percent in a range of W = !F c ±4t 
and E = 0±2t (d = 3, box distribution), while it acquires a strong energy dependence 
near the band edge. Generally, \u.(Dq/Dq) is a monotonically decreasing function as 
one moves from the band center to the band edge. The resulting energy dependent 
shift of (3(g) is again in agreement with exact numerical d iagonalizations [36]. A 







physical interpretation of the renormalization factor Dq/Dq would be required to 
make contact with the scaling behavior of the actual conductance g(L ) in the strong 
coupling regime, but has not been found yet. 

6. Conclusion 

A self-consistent theory of Anderson localization for the tight-binding model has 
been presented. The diagrammatic expansions in terms of Bloch and in terms 
of Wannier states were found to be equivalent, if maximally crossed diagrams 
are included self-consistently. Considering that no adjustable pararne- ters 
enter the theory, detailed quantitative agreement with results of exact numerical 
diagonalization of finite systems is found for all quantities available for comparison. 
The breakdown of the existence of a universal one-parameter scaling function at 
energies near the disordered band edge, as observed in the numerical calculations, 
can also be seen in our theory. 

The basic assumptions made in the diagrammatic theory presented here are (1) 
that at the Anderson transition there be no contributions which are not accessible 
by perturbation theory and (2) that there be no other infrared divergent diagram 
classes besides cooperons contributing to the p-h vertex. From the remarkable 
agreement with numerical data it appears, however, that, if such contributions exist 
- as suggested by field theoretical methods they should be important in a relatively 
narrow region around the critical point which is not resolved in our theory nor in 
the numerical calculations. This would require considering larger and larger system 
sizes. 
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